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^ Abstract 

t> 

Two-particle scattering in Ising field theory in a weak magnetic 
C\l field h is studied in the regime T > T c , using perturbation the- 

ory in h 2 . We calculate explicitly the cross-section of the process 
2 — > 3 to the order h? . To this order, 02^3 dominates citot (the 
probability of all inelastic processes) at all energies E. We show 
that at high energy the h 2 term in o" to t grows as I6G3 h 2 \ogE, 
rS where G3 is exactly the third moment of the Euclidean spin-spin 

correlation function. Going beyond the leading order, we argue 
that at small h 2 the probability of the 2 — > 2 process decays as 
E -idG 3 h 2 e ^00. 
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1 Introduction 



The Ising Field Theory (IFT) is the scaling limit of the 2D Ising model in a magnetic field H, taken 
near its critical point T = T c , H = 0. In field-theoretic terms, it can be defined as the unitary 
conformal field theory with c = 1/2, perturbed by its two relevant operators, 



Awt = A=i/2 cft ~ \ e{x)d 2 x + h I a(x)d 2 x . (1.1) 



2tt 



The field e(x), of conformal dimensions (1/2, 1/2), represents the temperature deviation from the 
critical point (i.e. m ~ (T — T c )); it is usually called the "energy density". The "spin density" a(x) 
has the dimensions (1/16, 1/16), and the associated coupling parameter h is the rescaled magnetic 
fielcQ Away from the critical point m = 0, h = the theory (1.1) is massiv^j and its physical 



content can be understood in terms of the spectrum of particles and their scattering amplitudes. 
The parameters m and h carry mass dimensions 1 and 15/8, respectively. Thus, apart from the 
overall mass scale, the theory (1.1) depends on a single dimensionless parameter r\ = m/|/i| 8 / 15 . 



The theory (1.1) admits an exact solution at two special ( "integrable" ) points. One is the case 



of zero magnetic field, h = 0, where (1.1) reduces to a theory of free Majorana fermions with 
mass \m\ (see e.g. fl|). This corresponds to the values r/ = ±oo, where the sign of the infinity 
distinguishes between two phases of the h = theory. At r] = +oo the spin-reversal symmetry is 
spontaneously broken, and the field a(x) develops a nonzero expectation value (cr(x)) = db a, with 
a known exactly [2j, 

a = s\m\ 1/8 , s = 2 1/12 e~ 1/s A 3 J 2 = 1.35783834. . . (1.2) 

{Aq = 1.28242712 ... is Glaisher's constant). At r/ = —oo the symmetry is unbroken, {cr(x)) = 0. 
Another solvable point is m = (i.e. r/ = 0), where the theory is not free, but an exact solution 
can be obtained due to the presence of infinitely many local integrals of motion [2]. 



At generic rj the theory (1.1) is not integrable. Nonetheless, its particle content is qualitatively 
understood since [3j. When r/ changes, the spectrum of stable particles gradually changes from 
an infinite tower of "mesons" in the "low-T" regime at rj — > +oo to a single particle in the"high- 
T" regime r/ — > —oo. On a quantitative level, various expansions of some of the masses near 
the integrable points are available through the perturbation theory [3l[7Tfl3], and rather accurate 



numerical data were obtained [7j|8j by means of the "truncated conformal space approach" of 14 
and by the numerical diagonalization of the lattice transfer-matrix |15[|16| . 

Apart from the integrable points, much less is known about the particle scattering. The numer- 
ical method of |14| is of limited power in addressing this problem. While this approach is successful 
in numerical evaluations of low-energy scattering phases [l7| , it is not clear how to apply it to an 
analysis of high-energy scattering. At the same time, understanding the scattering amplitudes, in 
particular their high-energy behavior, can shed some light on the structure of the theory. When 
the parameter rj changes, some particles lose their stability, becoming virtual or resonance states. 
Qualitative changes in the spectrum typically occur when the masses of some resonance states go 



1 We assume the standard CFT normalization of these fields: 

(e(x)e(0)) -> \x\- 2 , (<t(x)<t(0)} -> \x\~ 1/4 as \x\->Q. 

Also, note that the sign of m is changed as compared to the notations in [8][9j- 

2 Exception is the Yang-Lee point, which at real positive m appears at special purely imaginary h — 
± i m 15 / 8 (0.18933(2)). These are critical points [H], where the correlation length diverges, i.e. the mass of the 
lightest particle measured in the units of m vanishes. The associated CFT is c = —22/5 non- unitary Minimal 
model |6j. 
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to infinity (in the units of the stable particle masses). Analyzing this phenomenon clearly requires 
some knowledge about high-energy scattering. 

In this situation, perturbation theory around integrable points is a reasonable way to begin 
the analysis. In this work we study the 2 — > 2 scattering amplitude in the high-T domain m > 0, 
using perturbation theory in h to the leading order in h 2 . Unfortunately, in its present form, 
this perturbation theory is not as simple and transparent as the Feynman diagram technique. It 
relies on the intermediate-state decompositions, with the use of well-known exact matrix elements 
( "form- factors" ) of the field ct(x) at h = |18|. One of the technical difficulties of such "form- 
factor perturbation theory" is that separating disconnected contributions and the "external leg" 
mass corrections is not as straightforward as in the standard Feynman diagrammatic. We bypass 
these difficulties by using the optical theorem and the associated dispersion relation, expressing the 
amplitude through the total inelastic cross-section crtot- This procedure involves the well-known 
ambiguity in analytic terms, which can be fixed once the high-energy asymptotic of the amplitude 
is known. To eliminate the ambiguity we calculate this asymptotic directly (again, to the order 
h 2 ), using the techniques of It turns out that to this order the high energy behavior of the 
2 — > 2 S'-matrix element is dominated by the logarithmic term 

S 2 ^2 -> -l + 8/i 2 G 3 [log(E/E ) - iir/A] + 0{h 4 ) as E -> oo , (1.3) 

where E is the center of mass energy, and the coefficient G3 is exactly the third moment of the 
Euclidean spin-spin correlation function 

/•oo 

G 3 = r 3 {a(r)a(0)) dr = (2.3475292186... )m" 15/4 . (1.4) 
Jo 



The asymptotic behavior (1.3) indicates the logarithmic growth of the leading term in the h 



expansion of the total inelastic cross section. The unitarity bound otot ^ 1 suggests that (1.3) 
represents but the first terms of a series in the "leading logarithms" (h? \og(E)^ n . We present 
argument in support of this statement, the argument which also shows that the leading logarithms 

form the exponential series, resulting in the power-like decay | ^2 >2 1 ~ (E)~^ G ' ih2 at large E. This 

behavior suggests that while at low energies the two-particle scattering is dominated by the elastic 
2 — > 2 channel, at large E the scattering goes almost entirely into inelastic channels, even at small 
h 2 . 



2 General properties of 2 — > 2 S'-matrix element 

In this section we briefly describe the general properties of the 2—7-2 scattering matrix element, 
and the 2 — > n inelastic cross-section. For the most part, the content of this section is an adaptation 
of the textbook basics of the relativistic 5- matrix theory (see e.g. [19]), with minor simplifications 
specific for 1+1 kinematics. 

For simplicity, and in view of the problem at hand, we assume that the theory has one neutral 
particle, with mass m, which we refer to as particle A. The kinematic state of the asymptotic 
particle is characterized by its on-shell 2-momentum p^, conveniently parametrized by the rapidity 
9, p{9) = p^(9) = (mcosh#, msinh#). Here we use the notations 

\0l, 02, ■■■ , On)in(out) (2-1) 

for the -in (-out) state with n particles A, with rapidities 61, 62, • • • , n . We assume the following 
normalization of the one-particle states, 

(Q\0') = (2ir)5(8 -9') . (2.2) 



2 



In 2D, the 2—7-2 scattering is always purely elastic (i.e. the momenta of the two outgoing par- 
ticles are the same as the momenta of the incoming ones) by total energy-momentum conservation. 
Therefore one can write 



oo 



h, 02 )in — S2->2(Ql, 02)|fll; 02 ) out + ^ / ' ' ' x 



n=3 

\2 



X ^^d^iPin- P ou t)S2-,n(9u0 2 \(3 1 ,--- , /3 n ) \fa , • • • ,/3 n ) oui , (2.3) 

where P{ n = p(9i)+p(02) an d Pout = Y^i=xP{@i) are the total 2-momenta of the initial and the final 



states, respectively. Equation (2.3) defines the 2 — > n S-matrix elements S 2 ^ n - Our attention will 
be mostly on the element S^-^^i) $2)- Relativistic invariance demands that it actually depends 
on a single variable, the rapidity difference 9\ — 02] correspondingly, we will write 

S2^2(0i,02) = S{0 l -0 2 ). (2.4) 

The function S{0) has a direct physical interpretation at real values of (the physical domain of 
the s-channel scattering), but it can be analytically continued to the complex #-plane with certain 
singularities. The analyticity of S{0) was widely discussed in the specific context of factorizable 
scattering theory (see e.g. [20]), but much of this analysis goes through in the general case. The 
function S{0) is analytic in the strip < Qm9 < ir, except for possible poles at the corresponding 
segment of the imaginary-^ axis (such poles, if present, signify stable particles existing in the 
theory). S(9) takes real values at pure imaginary 9. The function admits analytic continuation to 
the whole 0-plane (with certain branch cuts, as is detailed below) via the functional relations 

S{0) = S(iir - 0) , (2.5) 

S(0)5(-0) = l. (2.6) 

The first of these relations expresses the crossing symmetry of the S'-matrix. The second follows 
from the unitarity of the 2—7-2 S'-matrix at the energies below the multi-particle thresholds. In 



the case of an integrable, purely elastic scattering theory, the multiparticle contributions to (2.3) 
are altogether absent, and as a result in that case S(0) is a meromorphic function on the whole 
complex plane. But in the general (non-integrable) case, of course, branching-point singularities 
associated with the multiparticle thresholds are present. There are branching points at = ±9^ n \ 
where 0^ are real positive solutions of the equation 

n , , 

cosh — = - ; (2.7) 



they correspond to the n-particle thresholds. By crossing symmetry (2.5), there are similar branch- 



ing points at the axis Qui 9 = -it, at 9 = in =p 9^; these represent the n-particle thresholds in the 



cross channel. Since the relations (2.5), (2.6) imply periodicity. 



S(6 + 2ivr) = S(0) , (2.8) 

the pattern is periodically extended along the imaginary-^ axis. We will define the principal sheet 
of the 0-surface by making branch cuts from flW to +00 and from —9^ to —00, with the periodic 
extension according to Eq.(2.8), as is shown in Fig. [T] (note the unusual directions of the real and 



imaginary axes in this Figure). Presumably, S(9) can be further analytically continued to under 
the branch cuts, to other sheets of the Riemann surface, but at the moment we do not have much to 
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Figure 1: Analytic structure of the two-particle scattering amplitude S(9) in the complex 0-plane. 
The solid lines represent the branch cuts associated with inelastic channels. The values of S(9) 
at different edges of the branch cuts represent physical 5- matrix element S, its complex conjugate 
S*, and the inverse values. The bullets • and circles o indicate possible positions of poles and 
zeroes, respectively. Poles located on the imaginary axis, within the physical strip < ^m9 < ir 
correspond to stable particles; poles on the strip — -k < Qm9 < are associated with resonance 
scattering states. 



say about its analytic structure there. Discontinuities of S(9) across the branch cuts are controlled 
by the probabilities of inelastic scattering events. As follows from the full unitarity condition, we 
have, at real positive 9 



S(9 + iO)S(-0 + iO) = S{9 + iO)/S(9 - iO) = 1 - a tot (9) 
where cr to t(0) is the total probability of all inelastic processes, 



(2.J 



O"tot(#) = ^ °2->n(0) , 



n=3 



(27T) 2 



sinh I 



~2n 



2tt 



(2.10) 



6 {2) (Pout " Pin) | S 2 -> n A,)f 



In the last line, the notations are the same as in Eq.( |2.3[ ) . Obviously, o"2-^n 

(0) = at \9\ < 0( n \ 
and hence the total cross-section (2.10) has the segment [6^ 3 \oo) as its support. 



Apart from these branch cuts, S(9) can have poles at the principal sheet. Due to the periodicity 
( 2.8 ) , one can limit attention to the strip — ir ^ 3m 9 ^ it. We will refer to the domain < 9 < ir 



as the "physical strip" (PS), since the values of S(9) at its boundaries represent physical scattering 



amplitudes. In view of equation (2.6), every pole of S(9) in the PS is accompanied by the associated 
zero in the strip — tt < Qm9 < 0, and vice versa, and for this reason we refer to the latter as the 



4 



"mirror strip" (MS). By causality, locations of possible poles within the PS are limited to the 
imaginary axis $te9 = 0. The poles in the PS are associated with the stable particles of the theory. 
In a unitary theory, a pole at = ia p (with < a p < tt), 



S(9) 



ir„ 



ia T , 



(2.11) 



with positive r p is the direct-channel manifestation of the stable particle ("bound state") A p with the 
mass M p = 2 cos(a p /2). Every such pole comes along with the cross-channel pole at ia p = i(TT — a p ), 
with the residue — ir p ^ On the other hand, poles in the MS are not restricted to lie on the imaginary 
0-axis. There may be poles at ia p 6 MS with real a p ("virtual states", in the terminology of the 
potential scattering theory), as well as with a complex a p (generally associated with resonance 
states). 

If the positions of all the poles on the principal sheet are known, the function S(9) can be 
written as 



s(9) = n 



sinh 9 + i sin a p 
sinh 9 — i sin a r , 



U(9), 



(2.12) 



where the product accounts for all the poles 9 = ia p , — a p ), including bound-state, resonance, 
and virtual ones. The function U(9) satisfies the same relations (2.5) and (2.9) as the amplitude 



S(9), but in addition it has neither poles nor zeroes on the whole principal sheet. It can be written 



as 



17(0) = exp{i sinh(0) A(0)} 



(2.13) 



with A(0) analytic everywhere in the PS. As follows from (|2.5|) and (2.6), A(0) satisfies the relations 



A(0) = A(ivr - i 
which in turn imply an enhanced periodicity 



A(0) = A(- 



A(0 + ivr) = A(0). 



(2.14) 



(2.15) 



Of course, A(0) has branching points at the thresholds 9^ n \ and in view of (2.9), its discontinuities 
across the associated branch cuts are related to the inelastic cross-sections. In particular, at positive 
real we have 

.log(l-a tot (0)) 



A(0 + «O) - A(0-iO) 



sinh 



(2.16) 



In view of Eq's (2.14), (2.15), it suffices to focus on the values of A(0) in the strip ^ 3rn0 ^ 
7r/2. The variable transformation 

u; = sinh 2 (2.17) 

maps this domain (with obvious identifications at the boundary) on the complex u>-plane with the 
branch cut along the real axis, from = sinh 2 0^ 3 ^ = 45/4 to +od^} We note that w relates to 
the center of mass energy E = 2mcosh(0/2) as 



w 



E 2 (E 2 -4m 2 ) 
Am 2 



(2.18) 



3 Of course, strict adherence to our assumption that there is only one kind of stable particle implies that there can 
be only two poles in the PS, one at ia.A = 2iri/3, and at i&A = in/3. 



4 It is useful to understand the mapping properties of (2.171. The edges Sm9 = +0 of the branch cuts 9 G 
[6^, +oo) and 6 € (— oo, — 8^] in Fig. Ill correspond to the upper and lower edges of the branch cut in the m-plane. 
The segment [0,6» (3) ] of the imaginary 0-axis is mapped onto the segment [0,-u/ 3 '] of the real u>-axis. Furthermore, 
segments [—1,0] and (— oo, — 1] of the real w-axis are the images of segments 5fte# = 0, QmO G [0,7r/2] and 9m9 = 
7r/2, 3?e# G [0, oo), respectively. 



5 



With some abuse of notations, let us write A(w) and ctot(w) for the quantities defined in Eq's (2.13) 



and (2.10), expressed in terms of the variable w. Equation (2.16) then relates the discontinuity 
of A(w) across the branch cut in the u;-plane to <Jtot(w), and one can write down the associated 
dispersion relation 

A(w) = A reg (w) + 1 / — — , (2.19) 

J w (3) {W-V)yjv 2-K 

where A reg (w) is an entire function of w, real at the real ■w-axia^J 

3 Scattering amplitude in IFT at weak coupling. 

Let us come back to specific theory ( |1.1[ ), and consider particle scattering in the high-T domain, 
at weak coupling \h\ <C m 15 / 8 . In this domain there is only one stable particle (3], which at h = 
becomes a free Majorana fermion, in particular S(0)\h=o = — 1q At small h the amplitude S(0) 
can be expanded in powers of h 2 . Here we mostly focus on the leading term of this expansion, the 
amplitude A(0) in 

To simplify notations, in what follows we set the units of mass so that 

m = l. (3.2) 

The amplitude A{0) is analytic in the PS, except for the poles at ictA = 2itt/3 and iola = *7r/3, 
representing particle A itself in the direct- and the cross- channels, respectively. Note that in view 



of the overall minus sign in (3.1), the residue of iA{9) at iaA must be positive. As follows from 



(2.5) and ( |2.6| ), A(0) satisfies the same symmetries as A(0) in (2.14), i.e. 

A{0) = A{-0) , A{0) = A(in - 9) , (3.3) 
and A(0 + in) = A(0). We will also see that 

A(0) = 0. (3.4) 



Again, writing A(w) for this amplitude expressed in terms of the variable (2.17), we conclude that 
this function is analytic in the whole to-plane, with the branch cut from uA 3 ** to +00, except for 
a single pole at w = — 3/4. The discontinuity of A(w) across the branch cut is related to the 
cross-section by the well known optical theorem (a simple consequence of (12.91) ) 



A(w + iO) - A{w - iO) = i^w <t£H(w) , (3.5) 

(2) 

where a\ Q [ is the leading term of the expansion 

a tot (w) = h 2 a£>(w) + 0(h 4 ) . (3.6) 



5 The integral can diverge if Otot approaches its unitarity bound 1 too fast at high energies; as usual, in such cases 
the appropriate subtractions are to be made. 

6 In 1 + 1 dimensions it is useful to distinguish between in- and out-states even in a free fermion theory. Thus, in 
theory (JTTTJ) with h = 0, the in-state with 61 > 62 > ■ ■ ■ > 6„ is identified with the Fock-space state \8i, . . . , 6 n ) = 
a'(0i)a' (62) ■ ■ ■ tr (0„)[O), while for the corresponding out-state the fermion creation operators must be arranged in 
the opposite order. 
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This leads to the representation 



w 



r w 



w + 3/4~^~ J w (3) (v — w) y/v 2ir 



dv 



(3.7) 



where r is a constant, yet to be determined. In writing (|3.7|) we have taken into account the 

echnically consistent 
oo; this asymptotic 



condition (3.4). We also did not add any analytic terms, which would be technically consistent 

y/w log it; as w 



with (3.5). We will show in Sect. 4 below that A(w) 



(together with (3.4)) rules out such terms. 



Of course, the amplitude A(0) admits standard perturbation theory representation, as the 
integral 

iA{0 l2 ) = -\ [ d 2 x(0 1 ,0 2 \Ta(x)a(O)\0 1 ,0 2 ) conn , (3.8) 



where integration is over Minkowski space-time, x = (x, t), and 0± 2 stands for the rapidity difference 
01 — 6 2 , which we assume to be non-negative. The integrand involves the connected parlQof the 
time-ordered product of the Heisenberg field operators. In principle, the matrix element in ( |3.8[ ) 
can be handled through the intermediate-state decomposition, with the use of the form-factors 



F(9 1 ,6 2 \p 1 , f3 2k+1 ) = (0i, 2 k(O)|/3i, • • • , fa+i) , 



(3.10) 



available in an explicit form |18|. Since the form-factors (3.10) vanish at 6\ — 9 2 = 0, this expan- 



sion directly confirms the property (3.4). However, this approach has well-known difficulties. In 



particular, individual 2k + 1-particle contributions are not Lorentz-invariant, this symmetry being 
restored only upon summing up all multi-particle terms. For this reason, it is easier to collect 



the multi-particle contributions by using the dispersion relation (3.7), expressing the amplitude 
through the h 2 term of the inelastic cross-section (3.6). The latter can be written in terms of the 



form- factors (3.10) as 



a 



E 

k=l 



(3.11) 



r( 2 ) 



(2fc + l)! 



/ II ^ (2vr)V 2 )(P in -P out ) \F{6 U 2 |/3i,..., (3 2k+l ) 



This leaves undetermined the coefficient r in the pole term in (3.7), which should be evaluated 
separately. 



3.1 Pole term 

The residue of A(6) at the pole = 2iir/3 can be deduced from the one-particle term in the 
intermediate-state decomposition of the matrix element in (3. 



, dPH^ W sinhff) ^ SMm)f + (312) 

cosn 0i + cosh 2 — cosh p + lU 



7 The connected part is obtained from the full matrix element by removing its disconnected parts 

(0i, 021^(20^(0)101,03)0,,™ = (0i,02|T(r(a;)a(O)|0 1 ,0i)-27r5(0i -0O(02|To-(2;)(t(O)|0 2 )- 

- 27r5(0 2 - 2 )(0i|T(r(a;)a(O)|0i) - (2tt) 2 8(9 1 - 9' 1 )S(e 2 - 02)(O!ro-(x)cr(O)iO) (3.9) 



before setting 9[ — 61 , 9' 2 = 02. In writing ( 3.9 1 we have omitted the terms involving <5(0i — 9' 2 ) and S(6i — 9' 2 ) (which 
are generally present by the antisymmetry of the matrix element with respect to 0i <-> 9 2 ) because they vanish at 

01 — 01) 02 — 02- 
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where . . . stands for the multi-particle contributions. With the explicit matrix elements from [18] 
we have 



(01, 2 1 <r (0)|/3) 1 2 = o 2 tanh 2 



coth 2 



coth 2 



(3.13) 



Since this expression is positive at all real values of the rapidities, we have erased the absolute 
value sign, making possible an analytic continuation to complex rapidities. The integral over /3 is 
eliminated by the delta- functiorj^J and the result explicitly shows a pole at = 9\ — 02 = 2iri/3, 



A(0) 



i 18\/3 a 2 
9 -2m/ 3 



(3.14) 



Note that the minus sign here appears because analytic continuation of the r.h.s of (3.13) to the 
pole point yields the real but negative value — tan 6 (27r/3) = —27. This, when combined with the 
overall minus sign in (3.1), Eq.(3.14) shows 2iri/3 to be a positive pole of 5(0), 

i36(ah) 2 



5(0) 



- 2vri/3 



+ 0(h q 



(3.15) 



The cross-channel pole at = in/ 3, with the opposite residue, can be extracted from the three- 
particle contribution to the intermediate-state decomposition, which also generates the threshold 
singularity at 0®. The higher multi-particle terms do not contribute to the residues, producing 
only the associated multi-particle threshold singularities. The residue in (3.14) corresponds to the 
value 



36 a* 



(3.16) 



of the coefficient in the pole term in (3.7). 



3.2 Three- and multi-particle contributions to a 



(2) 
tot 



to) 

The first term <r 2 4-3 m (3-11) can be evaluated in closed form (see Appendix A for details). When 
expressed in terms of the center of mass energy E = 2 cosh(0i2/2), it can be written as 

4%(E) = 0(E - 3) B(E) 1(E) , (3.17) 

where the step function Q(E — 3) is displayed in order to remind that this part of the cross-section 
vanishes below the three-particle threshold, 



B{E) 



4a 2 (E + 2) 



(2E-l) A (E-3f 



and 1(E) is the elliptic integral 

1(E) 
with 



(£-2)2 (E + 1)(E -1)2(E + 3)iE 3 



(3.18) 



1 



>t 2 



(i-xt 2 y. 



dt 



X 



(E + l)(E-3Y 
(£-!)(£ + 3) J 



(E-2)(2E + 1) 2 

^2 A > V 



(E + 2)(2E-l) 



E + 2 
E-2 



(3.19) 



(3.20) 



It is easiest to do this calculation in the center of mass frame 8i + 9-2 
is not Lorentz-invariant, the residue at the pole is. 



0. Although the one-particle term (3.121 
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Some basic properties of this expression are readily derived. Thus, near the threshold E = 3 

=.2 



a^ 3 (E) ^ Q(E - 3) 



5 6 Vf> a- 
2 5 3 4 V3 



(E-sy 



The integral (3.19) can be represented by the expansion 



which rapidly converges at any finite E. On the other hand, at E — > oo 



32a 2 r, ,„ 2 N \/3 7T-11 ,»v 
log(£ 2 ) + + 0(l/£ 2 ) 



7T 



(3.21) 



(3.22) 



(3.23) 



The behavior of a^ 3 (E) is shown in Fig. [2J 
40 - 




30 - 



20 - 



Figure 2: Energy dependence of the partial cross-section a^-ssiE) (solid line). The dashed line 



(2) 

shows the asymptotic form (3.23), which is seen to approximate 0"2_>3 very closely starting from 
relatively low energies. 

(2) 

The partial cross-sections a^^k+i with k > 1 are not evaluated in closed form. However, 



all values of energy E. Thus, the five-particle contribution a\ { 5 is smaller than 1% of <7 2 jU at all 



(2) 

numerical estimates show that these multi-particle terms in (3.11) are small compared to c^jU at 



,(2) 



values of E, and the higher terms in (3.11) are yet much smaller. Therefore, the amplitude A(8) in 



by a. 



(2) 

2—^3 V 



(3.1 ) can be approximated with high accuracy by the dispersion relation (3.7) with al Q l(v) replaced 



3.3 Amplitude 

(2) (2) 1 1 

With al ot (w) approximated by a^^iw), the second terms in (3.7) (which we denote A a (w)) can be 



evaluated numerically. The relative contributions of A a (w) and the pole term A p (w) = 36a w/(w + 
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Figure 3: Relative contributions of A p [w) (solid line) and A a (w) (dashed line) to the amplitude 
A(w) at low values of w. 



3/4) in (3.7) are shown in Fig. [3] for low values of w. One can see that in the shown domain 
of w the pole term brings a substantially larger contribution then A a (w). In particular, in the 
physical domain w > 0, but below the three-particle threshold = 11.25, A a (w) contributes 
no more then 15% to the scattering phase. For w > w^ 3 \ A a (w) develops the imaginary part 
yJw/A a^X^iw). At larger w the real part of A a increases, overtaking at w ~ 68 the pole term 
A p (w) in magnitude, and then approaches the asymptotic &a 2 yjw. At large negative w, A a (w) 
behaves as — (8<7 2 /7r) \J — w log(— w/wq), with wq = \ expjll — \/3 tt} rs 64.86. 



4 High energy limit of the amplitude 

In IFT with h = 0, the matrix element in ( |3.8[ ) admits exact representation in terms of the Lax 
equation associated with the special solution of the classical sinh- Gordon equation. In this section 
we describe this relation, and then use it to evaluate the asymptotic of ^(#12) at high energy, 

012 — > OO. 



4.1 The matrix element 



As is known since the classic paper [21], the two-point correlation functions G = (0\Ta(x)cr(0)\0) 
and G = (0\T fi(x) fi(0)\0) of the spin field a(x) and the disorder field in the h = IFT are 
expressed in terms of the special solution of the classical sinh-Gordon system 



d z dz¥ = - sinh(2v?) , d z dtx = q 
o o 

Here and below we use the light-cone coordinates 

z = x — t , z = x + t 



1 — cosh(2(/?) 



(4.1) 



(4.2) 



to label the points x of 1+1 Minkowski space-time. The relevant solution is Lorentz-invariant, i.e. 
the functions <p(z, z), x(z, z) depend only on the Lorentz-invariant combination p = zz = x — t 2 . In 
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this case (4.1) reduces to an ordinary differential equation, the famous Penleve III equation. The 
solution is specified by its p — > asymptotic form 



¥>0>) 



1 , P 
— log — 

2 S 4 



log(-O) + 0(p 2 Q 2 ) , X (p) = 7 log(16p) + log(-Q) + 0{p) , 



1 



where 



(4.3) 



(4.4) 



Under the special choice of the constant k 
p > 0, and decays as 



=7b , 



(which we assume) <p(p) is regular at all real 



<p{p) = -K Q (yfp) +0(exp(-3Vp)) 



7T 



(4.5) 



as p — >■ +oo [21] . The solution <p(p) can be regarded as the function of complex p. It is possible 
to show (using e.g. the large-/? expansions of Ref. [21]) that <p(p) is analytic in the whole p-plane 
with the branch cut from — oo to 0, except for the point p = where the logarithmic singularity 
(explicit in (4.3)) is located. Moreover, the asymptotic form (4.5) holds for — ir < arg(p) < tt as 
\p\ —7- oo. The function x(p) nas similar analyticity, except that at large \p\ it tends to a constant, 



X (p) = 4 log(s) + 0{e 2 ^~ p ) , |p| -»• oo , -vr < arg(p) < vr . 



(4.6) 



where s is the same as in (1.2). In the high-T regime ((a) = 0, (p) = a) of h = IFT the 



correlation functions G(p) and G(p) are written in terms of (p(p) and x(p) as follows, 

G(p) = e x ' 2 sinh (ip/2) , G(p) = e x / 2 cosh (<p/2) . 



(4.7) 



The correlation functions are real at space-like separations p > 0, while in the time- like domain the 
values at the upper edge of the branch cut, (i.e. G(p + iO) and G(p + i0)) have to be taken. 

This relation can be extended to the matrix elements of the products Ta(x)a(0) and Tp(x)p(0) 



sandwiched between any particle states, 
linear system, 

d 7 $>- = 



The main ingredient is a special solution ^±[x\P) of the 



and 



(4.8a) 
(4.8b) 

(4.9a) 
(4.9b) 



which constitutes the Lax representation of the sinh-Gordon equation (the first of Eqs.(4.1) guar- 



antees integrability of (4.8a) and (4.9a)), with e p playing the role of the spectral parameter. The 
solution ^±{x\f3) relevant to our problem is described in some detail in Ref. [9] [^J The components 
ty±(x\fi) are analytic functions on the double cover of Minkowski space-time parameterized by the 
coordinates x = (z, z), with branching singularities at the right and left parts of the light cone, i.e. 



9 In Ref. [9] this solution is described in Euclidean space, where z = x + iy, z = x — iy. Here we need its continuation 
y — > it to Minkowski space-time. The continuation is straightforward for space-like separations p > 0. In the time-like 
domain p < the standard iO-prescription z — > z + iOsign(z), z — > z + iOsign(z) is implied. 
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at z = and z = 0. The double cover is needed because *&±(x\/3) change sign when the point x is 
brought around the origin (z,z) = (0,0). 

Due to the obvious symmetry of Eqs.(4.8), (4.9), the functions i Sf±(x\(3) depend only on the 
combinations Z = e^z, Z = e~@z, and below we often use the notation ^>±(Z,Z) for them. The 
solution under consideration can be characterized by its (Z, Z) — )• behavior, 



vP+(Z, Z) 



7T 



-n \z 



Z\ 3 



7T 



-ft \z 



Z\ 3 



1 + 



1 + 



40 - 1 

40- 1 

~64~ 



Z l 



Z 1 



4 n 2\ 



64 



64 



Z l + 0(/> 4 



4r>2> 



Z l + 0(p 4 



(4.10a) 
(4.10b) 



where O is defined in (4.4). 

The matrix elements of the T-product a(x)a(0) between particle states are certain prod- 
ucts of the functions ^±(x\(3). The relation looks somewhat simpler for the "centered" product 
Ta(x/2)a(—x/2); this shift of course is irrelevant in the integral (3.8). First, define the two-particle 
matrix elements 



0(9,9') = (0\Ta(x/2)a(-x/2)\9,9'), 

G(9\e') = {9\Ta(x/2)a(-x/2)\9') - 2vr 5(9 -0')G, 



(4.11a) 
(4.11b) 



where G = G(p) is the two-point function (4.7), and we omit the argument x in Q(9, 9') and Q(8\6') 
to simplify notations. Then we have 



0(01,02) 
0(0l|0 2 ) 



G 



G 



+ e 
+ e 



* S (9 1 ,9 2 )-G* a (9 1 ,9 2 ) 



,02 



*a(01, 02) ~G 4^(01, 2 ) 



where G and G are the two-point functions ( |4.7[ ), and 

* S (01, 02) = Zi)tt_(Z 2 , 30 + Zi)tf + (Z 2j £ 2 ) 

*a(01, 02) = *+(^l, Z X )*_(Z 2 , ^2) " ^l)*+(^2, Z 2 ) 



Here 



(Zi,Z0 



(e Sl z, e" ei z) , (Z 2 ,Z 2 ) = (e 02 z, e^z 



)■ 



(4.12a) 
(4.12b) 



(4.13a) 
(4.13b) 

(4.14) 



Finally, the four-particle matrix element in (|3.8l) is the fermionic Wick product of the two-particle 



(4.15) 



matrix elements (4.12), 

(9 1 ,9 2 \Ta(x/2)a(-x/2)\9 1 ,9 2 ) i 



g- 1 [g(0 2 |0i)a(0i|0 2 ) - a(0i|0i)a(0 2 |0 2 ) + 2 )e(0i,0 2 )] • 



It will be convenient to write Q(0\9) (which is just the limit 9' = 9 of Q(9\9')) as the sum of two 
terms 

G(9\9) = G K(Z,Z)-G L(Z,Z), (4.16) 

where 



K(Z,Z) = * + (Z,Z)V_(Z, Z), 

l(z, z) = #_(z, z)d e ^> + (z, z) - * + (z, Z)^. (z, z) 



(4.17a) 
(4.17b) 



and again Z = e e z, Z = e 
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4.2 High energy limit of the integral (3.8) 



The expression for the matrix element in (3.8) described above is still too complicated for the 
integral to be evaluated in closed form. However, it is possible to use this representation to derive 
its limiting behavior at #12 — > 00. 

To simplify arguments, we assume that Q\ — > +00 and 2 — > — 00. In this limit the factors 
involving the ratio of e dl — e® 2 and e Sl + e 6 * 2 in (4.12) can be dropped, and the expression (4.15) 
simplifies as 

{e 1 ,9 2 \Ta(x/2)a{-x/2)\e 1 ,0 2 ) conn -> G \L{Z X , Z X )L(Z 2 , Z 2 ) + K{Z U Z 1 )K(Z 2 , Z 2 )\- 

- G \K{Z X , Z X )L{Z 2 , Z 2 ) + L(Z 1: Z{)K(Z 2 , Z 2 )], (4.18) 



where K and L are the functions (4.17). Furthermore, it is possible to show that at \Z + Z\ 1 
the functions *&±(Z,Z) approach plane waves. The exact form looks different for different regions 
in the (Z, Z)-plane. Thus, at Z — > +00, and Z>Zwe have 



#+(Z,Z) -> 2e v/2 cos 



Z -Z vr 

~~4 4 



(Z, Z) -> 2 e"^ 2 cos ( + - 

1 4 4 



(4.19) 



while at Z — > +00 and Z ^> Z 



*+(Z,Z) ^2e- Lp/2 cos 



Z - Z 



*_(Z,Z)^2e^ 2 COS { 2 ^ + * 

1 4 4 



(4.20) 



(when Z ~ Z — > +00 the difference between (4.19) and (4.20) is irrelevant since in this domain 
ip — > 0). Likewise, at Z — > —00 and 



* + (Z,Z)^2ie^ 2 cos(^ T ^ + ^ 



V-(Z,Z) -)• -2ie-^ /2 cos 



Z- Z vr 
4 4 



(4.21) 



while at Z - 

*+(Z, Z) 2i e^ /2 



00 and — Z ^> Z 

Z -Z 



cos 



7T 

+ 4 



*_(Z, Z) -> -2ie^ /2 cos 



Z- Z 



(4.22) 



In writing these asymptotics we have ignored the overall sign ambiguity of (if the sign at 
Z + Z -> +00 is fixed as in Eqs.( |4.19| ), ( |4.20[ ), the overall signs in ( |4.21[ ), ( |4.22[ ) still depend on the 
way of continuation around the point (Z, Z) = (0, 0)), since this ambiguity does not affect functions 



(4.17). From these equations one finds the corresponding asymptotics of functions (4.17), 

Z -2 s 



K(Z, Z)->2 cos 



L(Z,Z) ->\Z + Z\ 



(4.23) 



as \Z + Z\ — > 00. 

The leading 0\ 2 — > 00 asymptotic of the integral (3.8) clearly comes from the first term, 
G L(Zi, Z\)L(Z 2 , Z 2 ), in (4.18). The remaining terms involve the fast oscillating function K, and 
integration over (z,z) leads to their suppression by at least one power of e~® 12 , as compared to the 
leading term. Thus, the leading high energy asymptotic #12 — > 00 of the amplitude (3.8) is given 
by 



iA(9 



12) 



d 2 x\UV\G(p + iO), 



712 



+00 



(4.24) 
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where d 2 x = d~xdt = ^ dzdz, and 
U = Z 1 + Z 1 



e dl z + e~ ei z . 



V = z 2 + z 2 



e° 2 z + e'° 2 z . 



The integral in (4.24) can be handled as follows. Write 

\uv\ = uv-2uve(-uv) 



(4.25) 



(4.26) 



where 0(x) is the usual step function. The first term here is analytic in the coordinates (z,z), and 
its contribution can be evaluated using the Wick rotation t — > —iy, yielding 



d 2 xUVG(zz) — > 7rzG 3 exp(#i 2 ) as 6 , i 2 ->■ +oo 



where 



G, 



pG{p) dp 



(4.27) 



(4.28) 



is the third moment of the Euclidean spin-spin correlation function (1.4). The second term in 
(4.26) involves the step function which limits the integration domain to UV < 0. In the domain 
U < 0, V > (which lies entirely within the future light cone p = zz < , z < 0) one can 



parameterize z 



-p e 



-p e^, with the integration domain corresponding to the 



range 62 < 4> < 0\ oi the Lorentz parameter <fi. Therefore, at 612 — > +00 the contribution from this 
domain to (4.24) is 



UVG(p + iO) d 2 x 



1\2 exp(tfi 2 J 



t/<o<v 



pG(p + iO) dp 



612 exp(6>i 2 ) 



pG(p) dp, 



(4.29) 



where the last form is obtained by rotating the the integration contour p — > e~ l7T p, which is possible 
since G(p) decays exponentially into the upper half-plane of complex p. The "past" part of the 
light cone z > > z brings in an identical contribution. Collecting all these pieces together, we 
have 

j^|U-4G 3 [0u - iTr/2 + flb] , (4-30) 

where #0 is a rea l constant whose value can not be determined by the above simple analysis (its 
calculation would require much better understanding of the behavior of L(Z, Z) in the domain 
Z + Z-l). 



The asymptotic (4.30) is written in terms of the variable (2.17) as 



iA(w + iO) — > —2 G3 y/w [ log(w/wo) — iir] , w — > +00 



(4.31) 



where wq is a real positive constant related to 9q in (4.30). This equation can be extended to the 
whole complex w-plane with the branch cut along the positive part of the real axis, 

A(w) — > —2 G3 \/—w log(—w/wo), \w\ — > 00 , — 7r < arg(— w) < it , (4.32) 

where the branch of the multi- valued function is chosen in such a way that y/—w is positive and 
log(— w) is real at real negative w. From (3.5) one finds the high-energy behavior of the ~ h 2 term 
in the inelastic cross-section 



<4 2 J H 4G 3 log(«;/u;o) ~ 8G3 log(E 2 ) , 



(4.33) 
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(2) 

which is of the same logarithmic form as the high energy behavior of c^U' determined in sect. 3, 
Eq.(3.23). Note that the coefficient in (3.23) corresponds to the 1-particle contribution Gi-part (p) = 
(I/71") Ko(yfp) to the intermediate-state decomposition of the spin-spin correlation function G(p), 



pG , i_ part ( / o) dp 



Aa 2 



7T 



2.3475035314. . . 



(4.34) 



The smallness of the difference between the exact value of G3 in (1.4) and the one-particle contri 



(2) (2) 

bution reaffirms the statement that <r to j: is dominated by its three-particle component c^ 3 . 



5 Remark on higher orders in h 2 



The logarithmic growth (4.33) of the h 2 term in the inelastic cross-section, 

otot ^8G 3 h 2 log(£ 2 ) + 0(/i 4 ) as E -> 00 , 



(5.1) 



suggests that the higher-order terms of the h 2 expansion become significant at high energies, since 
unitarity requires Otot ^ 1 at all energies. One can expect that the terms of the order h 2n in the 
expansion (3.6) are ~ h 2n \og n {E 2 ) at large E. Here we argue that this is indeed the case, and 



moreover, that summing up these leading logarithms results in a power-like behavioi 10 

<7tot 1 - (E/E )- a 
I6G3/1 2 at small h 2 , and Eq 



(5.2) 

Const m, with yet unknown 



at E ^> 1, where the exponent a 
constant. 

First, let us note that the expression (4.24) admits a simple semiclassical interpretation, as 
follows. Let us introduce coordinates 



(u,v) = (e 6 * 1 z + e 



e 92 z + e~ 



-02 



(5.3) 



on the (z, z)-plane, and assume that at #12 —> 00 the particles with the rapidities 6\ and # 2 can be 
represented by their classical trajectories u = and v = 0, respectively. Also assume that unless 
the points x±, X2 are too close to the trajectories, the matrix element (61, Q 2 \Ta{x\)a{x2)\0\, 62} is 
approximated by the correlation function (0\Ta(xi)a(x2)\0), multiplied by the sign factor which 
depends on whether the insertion points x\ and X2 are located on the same side or on different 
sides of each of the trajectories. More specifically, we assume^ 



(Q\,Q 2 \To{x\)o{x 2 )\B\,Q 2 ) ~ sign(ui) sign(n 2 ) sign(u 1) sign(t> 2 ) {0\Ta(xi)a(x 2 )\0) , 



(5.4) 



where (u\,vi) and (^2,^2) are the coordinates ( |5.3| ) of the points x\ and x 2 , respectively^ 2 ! This 
expression still involves the "disconnected parts" [1 + (sign(«i) sign(vi) — 1) + (sign(u2) sign(w2) — 
1)] (0|Tcj(xi)(t(x2)|0). Subtracting those, one obtains 



46(-uivi)e(-«2«2) (0|T<r(a;i)<r(a;2)|0) , 



(5.5) 



{e 1 ,e 2 \Ta(xi)a(x 2 )\ei,e 2 ) CO v 

10 This behavior, with a = 64 (ah) 2 /(TY m 4 ), was earlier proposed by S. Rutkevich |23| on the basis of 
a different argument. He also conjectured the high-energy form of the partial cross-sections cT2->2+n ~ 
(l/nl) (a log(E/E )) n (E/E )-<*. 

xl To simplify notations, here and below we ignore the delta functions associated with the plane- wave normalization 
2.2|) of the particle states. Instead, one can th ink of normalized wave packets centered at the rapidities 81 and 82- 

in a different context. In our case it can be justified as follows. 



A similar approximation was used in Ref. 



22 



The correlation between the spin operators located at the points x% and X2 is due to exchanges by odd numbers 
of particles. An external particle trajectory which passes between these two points intersects an odd number of 
trajectories of these exchange particles, thus producing the minus sign. 
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where again 0(u) is the conventional step function. Since the correlation function (0\Ta(xi)a(x2)\0) 
depends only on the separation x\ — X2, one can explicitly integrate out one of the coordinates, 
1 

~2 . 

1 



d 2 Xl d 2 x 2 {9i,e 2 \Ta(x 1 )a(x 2 )\e 1 ,92}, 



2 sinh 9 
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d z x \UV\ (0\Ta(x/2)a{-x/2)\0) 



(5.6) 



where U = u\ — 112, V = v\ — V2 are the coordinates (5.3 ) associated with the separation x = x\ — X2- 



The last expression is equivalent to (4.24), which leads to logarithmic behavior (4.31). The origin 



of the logarithm is clear. Apart from the factor |Z7V|, the integrand in (5.6) is Lorentz-invariant 



The step function in (4.26) provides a cutoff for the integration over the Lorentz boost parameter 
<p in the integral over d 2 x 



dpdcj). The factor 9 



12 



^ log ti; in (4.29) is just the size of the Lorentz 



boost which brings the rest frame of particle 2 to the rest frame of particle 1. 

The same approximation can be applied to the higher-order terms in the usual perturbation 
theory expansion of 



S(9%2) = — {9\, 6*2|Texp { — ih / cr(x) d 2 x}\9 1 ,9 



2 /conn 



(5.7) 



by assuming that at #12 — > 00 one can write 

(9 1 ,92\Ta(xi)a(x 2 )...a(x2n)\9i,9 2 ) 



2/i 



Y[ sign(iij) sign(uj 



\Ta(xi)a(x 2 ) . . -Cr(x2n)\0) 



(5.8) 



where (ui,Vi) are the coordinates (5.3) of the points X{. 
sign(uj)sign(uj 



Again, apart from the sign factors 



the r.h.s. of (5.8) is invariant with respect to simultaneous Lorentz transforma- 
tions of the coordinates X\. Therefore, integration over n?=i d 2 Xi generates at least one factor #12. 
But it is easy to see that this integration in fact produces a much larger contribution ~ 0" 2 - Indeed, 
the 2n-point correlation function contains (2n)!/(2 n n!) fully disconnected terms 

(0|Tct(2; 1 )ct(x2)|0)(0|T ( j(x3)ct(x4)|0) . . . (0\T<r(x2n-i)<r(x2n)\0) + permutations. (5.9) 

These terms are invariant with respect to n copies of the Lorentz group, each acting on a respective 
pair of the points. Because of the factorization, these terms can be handled exactly as it was done 
above with the ~ h 2 contribution, and lead to the contributions ~ (h 2 #12)™ to (5.7). It is not 



difficult to see that the contributions from the above fully disconnected terms form exponential 
series in h 2 9\2 ~ h 2 logw. Assuming that there are no other sources of these "leading logarithms" 
(which of course needs to be proven by more detailed analysis), this suggests the power- like decay 
of the 2 



2 S-matrix element in terms of the variable (2.17), 

-2G3/1 2 



S(w) ~ (— w) ' ZU3tl " as \w\ -> 00 (-7T < arg(w) < tt) , (5.10) 
where again the branch of the function is such that S(w) is real at negative w. This corresponds 



to the behavior (5.2) of the inelastic cross-section. 



tion" of (3.1), 



Finally, we note that the asymptotic behavior (5.10) agrees with naive "two-particle unitariza- 
sinh6» + isin(2vr/3) sinh 9 - i sin(27r/3 + rh 2 ) ( ih 2 A a {9\ 



S(9) 



sinh# — zsin(27r/3) sinh# + isin(27r/3 + rh 2 ) 



exp 



sinh# 



(5.11) 



where r = 36a , and A a {9) is the second term in (3.7 
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13 Although (5.11 1 involves all powers of h 2 , it of course is not an exact expression for the 2 — > 2 ^-matrix element, 



since it ignores the higher-order corrections to the residue (3.151, as well as the higher order terms in <T to t 
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Appendix 

A Calculation of a. 



(2) 

To simplify notations, we will perform the calculations in the center of mass frame 9\ = —62 = 9/2, 
and assume that 6 > 0. It is useful to introduce notation x = e 9 ^ 2 and E = x + the last one 
has the meaning of the center of mass energy. Kinematics demands that E ^ 3, i.e. x ^ (3 + \/5)/2. 
We have 

x" 



w 



A 1 
2x 2 



The relevant 2 — > 3 matrix elements (3.10) have explicit form 18 

„2 x 



.r 



x 2 + 1 



G(x|{^}), 



with 



^/ 1 r in - A » + £fe 1 + xzk A z» - Zj 
G(x\{ Zi }) = a [ [ -^—y^ [[ ——f 

fc=l 1<J J 



(A.l) 



(A.2) 



(A.3) 



where z% = e^, i = 1, 2, 3. The expression (|3.11l) with = 1 takes the form 



a. 



(2) 



2 a 2 x 2 (x 2 - 1) 



where 
W{x) 



1 

3! 



2 ^ 3 vr (x 2 + l) 3 
cfei dz 3 



W(x) 



2a 2 VE 2 ^ 



7T 



E 3 



W(x) 



(A.4) 



*1 ^2 Z 3 



6{E - Zl -z 2 - z 3 )5{E - l/ Zl - l/z 2 - l/z 3 ) G z {x\{ Zi }) . (A.5) 



The integration in (A.5 ) is from to 00 for each of the variables zi. Since the integrand is obviously 
symmetric with respect to permutations of these variables, one can eliminate the factor 1/3! by 
limiting the integration to the domain < z\ < Z2 < z 3 < 00. The integral can be further simplified 
by changing to symmetric variables 



Sl = Z\ + Z 2 + Z 3 , S 2 = Z X Z2 + ZiZ 3 + z 2 z 3 



Z\Z 2 Z 3 



with the Jacobian 



D 



d(si,S2,s) 



d(zi,z 2 ,z 3 ) 



O2 - zi)(z 3 - z 2 ){z 3 - Zi) 



having the following form in terms of the variables (A.6), 

D = J s\sl - 4 s\s - 4 s\ + 18 sis 2 s - 27 s 2 



The integrations over s\ and S2 are eliminated by the delta-functions in (A.5), which set 

s\ = E , S2 = s E . 
At these values the Jacobian simplifies as 

D D(E\s) = 2El ^3(3 - s_)(s+ - s) , 



(A.6) 
(A.7) 

(A.8) 
(A.9) 
(A.10) 
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where s± are ordered (s + ^ S-) roots of the quadratic polynomial 



P(E\s) 



E A + 1SE 2 - 27 

IP 



8 + 1. 



(A.ll) 



Note that at E ^ 3 the roots are real and positive, and that they collide at the three-particle 
threshold E = 3. The remaining integral over s takes the form 



W 



s + G 2 (E\s) 
_ D(E\s) 



ds . 



(A.12) 



where G(E\s) is the factor G(:E|{,Zj}) in ( A. 2 ) expressed in terms of the variables ( A.6), and specified 

Iculatic 

G(E\s) 



to the values (A.9). Explicit calculation yields 

2E 2 + 1 R(E\s) D(E\s) 



E 2 -l Q(E\s) s 



where 



R(E\s) 



2 E(AE 2 - 7) 



Q(E\s) 



Es + 1. 



(A.13) 



(A.14) 



the integration variable to bring the limits to —1 and +1. We define the parameter 



To simplify analysis of the integral (A.12), it is useful to make a projective transformation of 

(A.15) 



A = X(E) 



1 



s+ + 1 



(E + l)(E-3)>- 



(E-1)(E + 3Y 



and trade s for the variable 



t 



1 1 



V\l + s 



This brings (A.12) to the form 

2(£ + 2) 2 (2£- l) 4 (£-3) 3 f 1 fl-fit 2 \ 2 V^t 2 



W 



(E-2) 2 (E 2 -1)(E + 3)I(E-1)I J-i \l-vt 2 ) (l-At 2 ) 



3-dt, 



(A.16) 



(A.17) 



where [i = n(E) and v = v(E) are given in Eq.( |3.20 ). Note that < fx < A < v < 1 at all E > 3, 
and at E — >■ oo these parameters approach 1. Combining (A.17) and (A. 4), one arrives at (3.17). 
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